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We report on results for the heavy flavor contributions to F2{x, Q^) in the limit Q^ 3> m^ at NNLO. By calculating 
the massive 3-loop operator matrix elements, we account for all but the power suppressed terms in m jQ . 
Recently, the calculation of hxed Mellin moments of all 3-loop massive operator matrix elements has been hnished. 
We present new all-A'^ results for the 0(n/ )-terms, thereby confirming the corresponding parts of the 3-loop 
anomalous dimensions. Additionally, we report on first genuine 3-loop results of the ladder-type diagrams for 
general values of the Mellin variable TV. 



1. Introduction 

Deep-inelastic scattering (DIS) in the region of 
large enough values of the gauge boson virtu- 
ality Q^ — ~q^, allows to measure the leading 
twist parton densities of the nucleon, the QCD- 
scale Aqcd, and, related to this, the strong cou- 
pling constant as{Q^) = as(Q^)/(47r), to high 
precision, [I] . Especially in the region of smaller 
values of Bjorken-x, the DIS-structure functions 
F2,l{x,Q^) contain large cc-contributions of up 
to 20-40 %, therefore deserving further investi- 
gation. Our goal is the calculation of the com- 
pletely inclusive heavy flavor Wilson coefhcients 
which constitute the perturbative contributions 
to the structure functions and are denoted by 
Hj,(2,L)(a;, Q^/At^,m^//x^). Here, we include ef- 
fects of one species of heavy flavors in the final 
state and/or its virtual contributions. At present, 
the heavy flavor Wilson coefficients are known ex- 
actly to NLO by a semi-analytically result in x- 
space [2] . Due to the size of the heavy fiavor cor- 
rections, it is necessary to extend the description 



of these contributions to 0{al), and thus to the 
same level which has been reached for the mass- 
less Wilson coefficients ^. 

A calculation of these quantities in the whole 
kinematic range at NNLO seems to be out of reach 
at present. However, in the limit of large virtu- 
alities, Q^ > 10 m^ in the case of -F|'^(x, Q^), one 
observes that F^^ix^ Q^) are very well described 
by their asymptotic expressions '4' disregarding 
power corrections in ni? /Q^. In this kinematic 
range, the heavy flavor Wilson coefficients have 
been obtained analytically for ^^'^(x, Q^) to 2- 
loop order in |4l5j and for Ff^lx.Q"^) to 3-loop 
order in [6]. The asymptotic expressions are ob- 
tained by a factorization of the heavy quark Wil- 
son coefficients into a Mellin convolution of mas- 
sive operator matrix elements (OMEs) Ajk and 
the massless Wilson coefficients Cj^i. In case of 
only one heavy flavor species this factorization 
reads in Mellin-space, |4|7j , 



* Speaker 



H„(2,L)(iV) = A,,(7V)-C,,(2.L)(A^), 

i,j = q,g ■ (1) 

Here, we indicated the dependence on Mellin-A^ 



and suppressed all further variables, of. [418] for 
details. Eq. ([T]) allows to calculate the heavy fla- 
vor Wilson coefficients in the limit Q^ ^ m? up to 
0{al) by combining the results obtained in Ref. 
[3] for the light flavor Wilson coefficients with the 
3-loop massive OMEs. 

Another physical motivation for this calculation 
is that the massive OMEs also serve as transition 
functions if one wants to define parton densities 
for massive quarks in the framework of a variable 
flavor number scheme, 0- This is of particular 
interest for heavy quark induced processes at the 
LHC, such as cs -^ W~^ at large enough scales 
Q^ . Another important point is that in the course 
of our calculation, we obtain the 0(n/ )-parts of 
the 3-loop anomalous dimensions, [9^, which are 
thus confirmed for the first time in an indepen- 
dent calculation. Finally we are also interested 
in the mathematical structures of the Feynman- 
parameter integrals emerging in our calculation, 
leading to new insight, cf. [lO] . 
In the following, we briefly describe the calcula- 
tion of the flxed moments of the massive OMEs 
and then present our recently obtained all-A^ re- 
sults. For more details on the challenges we en- 
countered on the mathematical side, see |10j . 

2. Massive OMEs and Fixed Moments 

We are interested in the flavor-decomposed 
twist-2 massive OMEs 
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The external on-shell particles are denoted by « = 
q, g and Ofc stands for the quarkonic (fc — q) or 
gluonic (fc = g) composite operator emerging in 
the light-cone expansion. The subscript H indi- 
cates that we require the presence of heavy quarks 
of one type with mass m. fi'^ denotes both the fac- 
torization and renormalization scales. The com- 
posite operators give rise to additional Feynman- 
rules depending on Mellin-iV which can be found 
in Ref. 0. 

In case of the gluon operator, the contributing 
terms are denoted by Agq^Q and Agg,Q. For the 



quark operator, one distinguishes whether the 
current couples to a heavy or light quark. In 
the non-singlet (NS)-case, the operator, by def- 
inition, couples to the light quark. Thus there 
is only one term, A^q n- In the singlet (S) 
and pure-singlet (PS)-case, two OMEs can be 
distinguished, {AHq, A^gQ} and {A^(fq, A%g}, 
where, in the former case, the operator couples 
to a light quark and in the latter case to a heavy 
quark. We also consider the transversity opera- 
tor, giving rise to the OME AtA^|q . 
Up to and including 2-loop order, all massive 
OMEs have been calculated in Refs. [4|7j and 
conflrmed independently in j5|ll| . The transver- 
sity terms have been obtained in Ref. [12]. Ad- 
ditionally, all 2-loop 0(e)-terms in D = 4 -|- e 
dimensions have been calculated in |11|12|13] . 
These terms are needed for the renormalization 
of the massive OMEs at 3-loops. Let us briefly 
review the most important steps of renormaliza- 
tion, j4 j8l . We apply the MS-scheme, except for 
mass renormalization, which is performed in the 
on-shell scheme [T^. After mass and coupling 
constant renormalization, the remaining singu- 
larities are of the ultraviolet and collinear type. 
The former are renormalized via the operator Z- 
factors, whereas the latter are removed via mass 
factorization and absorbed into the parton densi- 
ties. Note that in the last two steps the anoma- 
lous dimensions of the twist-2 operators emerge. 
Thus at NNLO the fermionic parts of the 3-loop 
anomalous dimensions calculated in Refs. |9I15) . 
cf. also [16], appear. The general structure of 
the un-renormalized and renormalized OMEs at 
3-loops is 
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The pole terms in e of Eq. ([3]) and the logarithmic 
terms in vr? j ^ in Eq. ([4]), respectively, are com- 
pletely determined by renormalization and can be 
expressed in terms of the anomalous dimensions 
up to 3 loops, the expansion coefficients of the 
QCD /3-function up to 2 loops and the 1- and 2- 



loop contributions to the massive OMEs. For the 
exact formulae, we refer to ^8 . Hence these terms 
are already known for general values of N and can 
be used for first phenomenological analyses, |17j . 
This is not the case for the constant term, which 
contains the genuine 3-loop contributions a^-- . 
The massive OMEs at 0{al) are given by 3-loop 
self-energy type diagrams, which contain a local 
operator insertion. The external massless parti- 
cles are on-shell. The heavy quark mass sets the 
scale and the spin of the local operator is given 
by the Mellin-variable N. Using the programs 
QGRAF, [18 , color, [TH], MATAD, [10], and the 
computer algebra system FORM, [31], we calcu- 
lated fixed moments of all OMEs. For the terms 
Qa' qg Q ^^'^ gg Q *^^ even moments N ^ 2 
to 10, for AgJ'^^ to iV = 12, and for ^^Jq^' 
^iJ'n I ^nJn to N = 14 were computed. For 
the flavor non-singlet terms, we calculated as well 
the odd moments iV = 1 to 13, corresponding to 
the light flavor (— )-combinations, 0- We also 
calculated the transversity OME IS.tA 'q for 
the even/odd moments N = 1..13, [12]. All our 
results agree with the predictions obtained from 
renormalization, providing us with a strong check 
on our calculation. For all moments we calcu- 
lated, we agree with the known results in the lit- 
erature, especially with the T/~parts of the 3- 
loop anomalous dimensions, [9115116] . Our com- 
plete results are given in Refs. [8112] and will 
be used for first phenomenological parameteriza- 
tions, [T7] . 

3. The n/ contributions for all N 

As a first step towards the calculation of the com- 
plete 3-loop massive OMEs for general values of 
N , we consider the ny-contributions only. In the 
following we describe the computation of these 
terms and present our results for the constant 
term aQ and the 3-loop anomalous dimension 

7qg which enters the single pole term in Eq. (jSj). 
The remaining OMEs will be given in [22j , cf. also 



In Refs. |5lllll3j . the 2-loop massive OMEs were 
calculated using Feynman-parameters in order to 
arrive at a representation of the momentum inte- 



grals in terms of finite and infinite sums which 
depend on the Mellin-variable N . The summand 
is a hypergeometric expression in terms of F- 
functions, which could safely be expanded in e, 
before performing the summation. Thus we ob- 
tained up to twofold sums, one of which derived 
from a hypergeometric function 3_F2 extending to 
infinity and the other usually being a finite sum. 
These sums could then be calculated applying al- 
gebraic or analytic summation techniques and for 
the more complicated sums the summation pack- 
age Sigma, [10|13|24j . 

For the 3-loop n/-terms, we followed in princi- 
ple the same approach. In this case an additional 
light quark loop is present along with the heavy 
quark loop, compared to the 2-loop case. There 
is one basic topology only, which is shown in Fig- 
ure [TJa). From this all contributing diagrams 
can be derived by attaching outer light partons 
or ghosts and inserting the quarkonic operator in 
all possible ways. Since the external particles are 
on-shell, for the lowest moment all diagrams can 
be reduced to this massive tadpole. By recur- 
sion, this holds for higher moments as well. As 
an example, in Figure[Tl (b), one of the more com- 
plicated diagrams contributing to the nj-part of 
A^qI is shown, cf. [10]. 
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Figure 1. Examples for 3-loop diagrams oc n/. 
Solid lines: heavy quarks, dashed lines: light quarks, 
curly lines: gluons. ®: operator insertion. 

The calculation of the massive tadpole shown in 
Figured! (a), is straightforward. After integrat- 
ing out the light quark loop, one is left with a 
2-loop massive tadpole, for which the result is 
known analytically for arbitrary exponents of the 
propagators. However, as in the 2-loop case, the 
calculation is significantly more complicated once 
the operator insertion is present. Nevertheless, all 



momentuni integrals could be re-written in terms 
of up to two fourfold sums. These sums differed 
in their structure from those encountered in the 
2-loop case, making the calculation significantly 
more complicated. Here, we refer to [iQi for more 
details on how these sums were calculated using 
the MATHEMATICA-based program Sigma, [24] . 
Let us now present some of our results. For the 
""^-term of the 3-loop anomalous dimension ^qj 



we obtairo 
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with the polynomials 

Pi = 57V'' + 2QN^ + 47A^2 + 437V + 20 , (6) 
P2 = ^N^ + 20Ar^ + 41A^2 _^ 49^ _|. 20 , (7) 
Pz = 197V^ + 124iV^ + 4927^ + 1153iV3 



+1362iV^ + 712iV+152 , 
P4 = 1152 + 7296iV + 19904Ar2 



(8) 



^Note that the term oc n'J in _9_ corresponds for the mas- 



sive OMEs to the term oc n 



+30864iV3 + 258967^ + 68007V5 
-10364iV^ - 8899iV^ + 3557iV^ + 8534Ar^ 
+5194iVi° + 14857V" + IQbN^'^ , (9) 

P5 = 47Ar^ + 145^^'"^ + 4267V2 + 412iV 
+120 , (10) 

Pg = 997V" + 9907V" + 49257Vi2 
+179167V" + 466497V^° + 724467V^ 
+322837V* - 955927V^ - 2675247V^ 
-4794727V5 - 5869287V^ - 4551687V3 
-269760/V2 - 1221127V - 27648 . (11) 



The above and the following expressions are given 
in terms of harmonic sums, [25126) . taken at argu- 
ment TV. Note that we perform an analytic con- 
tinuation to complex values of 7V starting from 
the even moments. Eq. ([5]) agrees with the cor- 
responding results in the literature for the case 
of general values of 7V and for fixed moments 
[8I9I16| . Our new result is the constant term, 
which reads 
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The terms Qi . . . Q20 are rather lengthy irre- 
ducible polynomials in N which will be given in 
Ref. [22]. Eq. p^ agrees for fixed moments with 
the results of Ref. [8]. The term Aq forms the 
most complex contribution to the 0{nf) terms of 

(3) 

qg,Q' 

PS, (3) 



the massive OMEs, along with the term A 
The latter as well as the 0(n/)-terms of A 



.PS,(3) .NS,(3) 



and At^„„ 'A have been ob- 



tained by us as well and will be presented in (55] 
in detail. Furthermore, the corresponding con- 
tributions to the remaining 3-loop anomalous di- 
mensions in the vector- and transversity case are 
obtained. In this way, an independent recalcula- 
tion of these quantities given in |9I15I16) before, 
was performed using very different methods. 
Finally, it is interesting to consider the small- 
and large-TV limits. We obtain for the n^-terms 



of the renormalized OME Aq 
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ln{N) = ln{N) + 7^ , (15) 

with -fE being the Euler-Mascheroni constant. 
The ^/-contributions at 0{a^) to the massive 
OMEs contain nested harmonic sums up to 
weight w = 4. This also applies to all individ- 
ual Feynman diagrams, which we calculated in 
Feynman-gauge. After reducing the contributing 
harmonic sums to their common basis using alge- 
braic relations [27^ between them, the following 
harmonic sums contribute 

Si, S2, S-2, S3, S-3, 

^2,1, 5-2, li S4, 5-4, 

53,1, 5-3^1, 5_2,2, "5*2,1,1, 5-2,1,1 • (16) 

Note that this set of harmonic sums does not 
contain the index { — 1}. Due to the structural 
relations [28130) . differentiation and argument- 
duphcation, cf. [26], the set ([T6| can be reduced 
even further. 5i represents the class of all single 
harmonic sums. Hence only the six basic har- 
monic sums 

5i ,52,1, 5'_2,i ,5'_3^i, 5*2,1,1, 5'_2,i,i . (17) 

are needed to represent the 3-loop results for the 
ny-contributions to the OMEs. 
In intermediate steps, we also observed so-called 
generalized harmonic sums [29130] . They obey 
the following recursive definition [29] : 

5,„i^...(a;i,...; A^) = 

/ , Trr[ / J ■m2 '-'"13,... [X3, ..., *2J 



11=2 



i2 = l 



f*^r?ii+m2,Tn3,... 1*^1 ' '^2j -^Sj ---^ -^^ ) 



(18) 



In the present calculation the values of Xi extend 
to {—1/2,1/2,-2,2}. These sums occur in lad- 
der like structures, cf. [3122131] and the next sec- 
tion. They may also emerge, however, if contri- 
butions to 3-loop Feynman diagrams containing a 
2-point insertion are separated into various terms 
in an arbitrary way. The weight of these sums can 
reach w = 5, although only w = 4 sums should 
remain in the final results. Examples for these 
sums are : 

§i(l/2;iV), S2{-2-N), 524(-l,2;iV), 
§3,i(-2, -1/2; N), ^2,3(-2, -1/2; N), 
^i,i,i,2(-l,l/2,2,-l;7V), etc. (19) 

The algebraic and structural relations for these 
sums are worked out in Ref. 30i. Similar to the 
case of harmonic sums, corresponding basis rep- 
resentations are obtained. Expressing our results 
in terms of this common basis, all generalized 
harmonic sums canceled for each individual di- 
agram. In the beginning, however, it has been 
unclear, whether this was bound to happen. It 
is not excluded that these sums contribute in the 
final results of some of the massive OMEs at 3- 
loop order. 

4. Ladder type structures 

As the next step towards a complete calculation of 
the 3-loop OMEs, we consider diagrams deriving 
from the ladder-type topology shown in Figure 
121 cf. [H]. Contrary to the diagram shown in 
Figure [Tfa), it can not trivially be reduced to a 
2-loop diagram. 
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xlx"^ = I dxi 
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(1 - 2:1X3 - X2Xi)"- 



(20) 



for arbitrary exponents of the propagators. This 



leads to a double infinite sum (m; 



'Vj^ etc. 



h oc 



y^ (-2 - g/2 + t^l2),n(-2 - g/2 + Vi5)n 

(2 + £/2)„,+„(2 + m + e/2)^,, {2 + n + g/2)_,, 
m\n\ 

(21) 

with {a)c = r(a -|- c)/T[a) the Pochhammer sym- 
bol. We checked this representation for various 
integer values of the Vi using MATAD and found 
complete agreement. Note that for all diagrams 
deriving from this topology, the Fi -structure oc- 
curs due to the diagram's topology and mass dis- 
tribution, and its form is independent of the op- 
erator insertion. As in the previous section and in 
the 2-loop case, this allows to obtain a represen- 
tation in terms of a multiply nested sum, which 
is regularized and can be expanded in g, for each 
diagram belonging to this class. Consider as an 
example the scalar diagram shown in Figure [S] 
with all Vi — 1. 




Figure 3. Example 3-loop ladder diagram. 



Figure 2. 3-loop ladder tadpole, 
nents of the propagators. 



denote expo- 



However, one finds that in the scalar case it can 
be represented in terms of an Appell-function of 



One obtains the following parameter integral 

h = Cz I dxi xl (1 - X3)2a;| 
Jo 

(1 - X4)4(x5(l - x^ +a;4a;6(l - Xi - x^) 

^Xi^XxX'j -I- X^X2X^) ^ X]^"^ X2^ 6'(1 — X\— X2) 



(1 -Xi -X2) 



a;3 - 1 X4 - 1 

1 - Xi X2 



X3 



X4 



-2+3e/2 



(22) 



C3 is a trivial pre-factor. The operator insertion 
contributes as a polynomial linear in each Feyn- 
man parameter raised to the power N. This term 
is in a certain sense "convoluted" with a structure 
deriving from the Fi in Eq. (|20|) . The application 
of the binomial theorem implies additional sums, 
until the integrand isgiven in such a way that Eq. 
(f20| can be applied Q Thus one obtains: 



h = 
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[N + l){N + 2) 



00 00 N + 2 /,T . r,\ I 
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j 



J=2 



U'=l ^ ^ r=Q ^ ' 



(fc, 
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{j + r + I + m + n - e){N + 3 - j) 



B(j,m + 1 + - 



(23) 



Here, i?(a, b) is the Euler-Beta function. After 
expanding in e, this sum can be reduced by Sigma, 
cf. [TOl and references therein. One obtains 



h = 
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{N + 1){N + 2){N + i) 6 




7V2 + i2Ar+l6 



Si 



4(2iV + 3) 



2(iV+l)(iV + 2) ' {N + IY{N + 2) 



Si 



^At first sight, the last line of Eq. H22| l and the denom- 
inator of Eq. I I2OI I seem to be different. After applying 
analytic continuation relations for Fi, their structure be- 
comes the same 
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-TT^l 
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^3-2(-l)^'5-2,i-(iV + 3)52, 
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2™+*5'i,2 7:,1;A^ 



-2^+'^i,i.i(^,l,l;iV) }+0{e) 



(24) 



It is interesting to note, that in this case gen- 
eralized harmonic sums appear in the result for 
one scalar integral. For the n/-terms, this never 
happened. There, the generalized sums emerged 
only when splitting up the complete contribu- 
tion to one diagram into several pieces, without 
further reference to certain momentum integrals 
anymore. This could point towards generalized 
harmonic sums remaining for individual diagrams 
or even for the complete physical OME. Another 
interesting point is that the powers 2^ lead to di- 
vergences as N —^ 00 and are therefore expected 
to cancel in the physical expression. Details on 
how to calculate this class of ladder-type dia- 
grams and results for the most complicated scalar 
integrals will be given in [31j . 



5. Conclusions 

The heavy flavor 3-loop Wilson coefficients are 
needed in consistent analyses of the DIS World 
data at NNLO. In case of the structure function 
F2{x,Q^), the asymptotic representation applies 
for scales Q"^ jw? > 10. After a larger number of 
Mellin moments were computed, now the asymp- 
totic massive Wilson coefficients are calculated 
for general values of the Mellin variable N. Here, 
we studied both the bubble- and ladder topolo- 
gies and developed the corresponding computa- 
tional frame. A first class of contributions has 
been completed in the case of the quarkonic oper- 
ator insertions with all contributions to the color 
factors TpnfCA,F- 
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